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Abstract. This paper aims to demonstrate the applicability of the £2-integral 
transform to Partial Differential Equations (PDEs). Of special interest is sec- 
tion (6), which contains an application of the £2-transform to a PDE of expo- 
nential squared order, but not of exponential order. Sections (1) and (2) aim 
to introduce the history and some elementary properties of the £2 transform, 
(3) and (4) include some of the transform's simple applications, and section 
(5) introduces the £2 convolution. 



1. Introduction 

Integral transforms, such as the Fourier and Laplace transforms, have numer- 
ous applications to differential and partial-differential equations. For example, the 
Fourier integral transform can be used to solve the popular heat equation, vari- 
ations of which form the basis for the partial differential equation that leads to 
the Black-Scholes model for options pricing in Europe. Alternatively, the Laplace 
transform has applications to physics in the analysis of linear-time invariant sys- 
tems, like electrical circuits and mechanical systems. The /^-transform already has 
some well understood applications to classical ordinary differential equations such 
as Bessel's differential equation [2](Yurekli, Wilson (2002)) and Hermite's differen- 
tial equation[3j(Yurekli, Wilson (2003)). The goal of this article is to illustrate the 
applicability of the /^-transform to solving certain partial differential equations-an 
application of the /^-transform never before examined in literature. Of particular 
interest in this paper is an example of a partial differential equation whose solution 
can be found by using the ^-transform, but is not of exponential order. This fact 
is significant in the sense that the nature of the solution assures us that it could not 
have been achieved through use of the Fourier transform, which can only handle 
equations of up to exponential order. 

2. Properties and Definitions Related to the £ 2 Transform 
Recall that a function / : [0, 00) — > R is called exponential squared order if 

lim f(x)e~ x2 = 

x— >oo 

Definition 1. For any exponential squared order function fit), the ^-transform 
of f is defined as: 

oo 

2 „2 



C 2 {f{x)-s}= / xe- x s f(x)dx 



o 



Date: February 16, 2012. 



1 



TODD GAUGLER 



Example 1. For n > the following is true: 



(1) C 2 {x 2n ;s} = 



2 s 2n+2 

In the case where n = 0, the definition of the C 2 -transform gives us the following: 

„2„2 1 

2? 

When n — 1, we arrive at the following through integration by parts: 



C 2 {l;s} = / xe x s dx - o , 



Property ([1} then follows by a simple induction. 



There exists a differential operator 5 X , denned as follows: 

S x = -- — 

x dx 

We now recall some of the properties of the /^-transform. ( pQ, [2]). 
Proposition 1. Let f be a function of exponential squared order. For all n > 0, 

(2) C 2 {S x f(x); s} = 2s 2 C 2 {f{x)- s} - /(0+) 
where /(0 + ) = linij.+^o f{ x ), an d 

(3) C 2 {x 2n f(x); s} = ^2-6?£ 2 {f(x); s}. 
Proof. For the first claim, we calculate 

poo i J poo 

xe' x2s2 8 x f(x)dx = / xe~ x2s2 ■ - — f(x)(dx) = / e^' 2 f {x)dx. 
to Jo x dx J 

Integrating by parts, we get 

f(x)(e- x2s2 ) + / 2xs 2 xe- x2s2 f(x)dx. 
o Jo 

Evaluating from to oo and using the fact that / is of exponential squared order, 
we can write this expression in terms of the /^-transform: 

2 S 2 £ 2 {/(x); S }-/(0+). 

For property © , taking the case in which n=l, we get: 



S s C 2 {f(x); s}=6 s xe~ x s ~ f(x)dx 
Jo 

We can bring the differential operator 5 S inside this integral, we obtain: 

d 2 2 f°° 2 2 

— • xe~ x2s2 f(x)dx — -2 / x 2 ■ xe~ x2s2 f(x)dx 
ds J 

We remark that the last term of this equation is equal to: 

£ 2 {-2x 2 f(x);s} 

For n > 1, through induction we arrive at 

5 n s C 2 {f{x)-s} = -2 n C 2 {x 2n f{x) ]S } 
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Written alternatively, 

On s 



C-l ) n 

£ 2 {x 2n f(x);s} = [ —J-S:C 2 {f(x); S }. 



a 



3. A First Application 
Consider the following partial differential equation : 

(4) t 3 u tx + 2xu = and u(0 + ,i)=0 

where it = i) for a;,t > and u(0 + ,t) = lim u(x,t). 

Writing the PDE in Equation §4§ [ n terms of the differential operator, we have: 

,1 d 

t s — —u t + 2u = 

x ax 

or equivalently, 

t 3 5 x ut — —2u. 

Taking the /^-transform of both sides and using property ([3]) we obtain 

2s 2 t 3 u t - u(0 + ,t) = -2u 
where u — u(s, t) = C 2 {u(x, t); s}. Equivalently, 

-1 „ u(0+,t) 

By the initial condition in (QJ, this last term is zero. One solution to this differential 
equation is 

1 

2? 

We now write u(s, t) as a series to obtain: 



1 oo f t~ s~ \n 00 1 ? 



2s 2 ^ n\ ^ 2s 2n + 2 ■ t 2n ■ n! ' 

n=0 n=0 



Using property ([T]), we can calculate that u(x,t) — C 2 1 {(u(s, t); x} is 



t 2n (n\) 2 ' 

n=0 y 1 

4. Further Generalizations 
Similarly, we can solve all partial differential equations of the form: 

(5) = f(t)u + f{t)-u x + g{t)-u xt 

x x 

Taking the C 2 -transform of the PDE in ([5]), we compute: 
= f(t)u + f(t)2s 2 u + g{t)2s 2 u t 

Rearranging: 

itt (l + 2s 2 ) 



u 2s 2 



M(t) 
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Where M(t) = 4$. Claiming that L(s) is any function of s, a solution for u is as 
follows: 

u( s , t) = e^^'/o M ^ . L(s) = do . e ^P^ . L {s) 

Writing this as a series-denoting J* M(w)dw as M(t), and letting L(s) = ^j- we 
get: 

-M(t) n ^ 7W(t) n ^ M(t) n ^ M(i) n 



-s 2 ^ n! • (2s 2 )™ ^ ra! ~ ^ n! • (2™s"+ 2 ) ' ^ 

n=0 v 1 n=0 n=0 y 1 n=0 

M(t) n rS, 



n i . 2(«-i) . (2s n+2> ) ^ 

n=Q ' v 7 



n=0 

using identity @, and taking the inverse ^-transform, we claim that: 

v ' ; ^ n!) 2 • 2"- 1 ^ 2n! 
Lastly, we can always solve the following types of partial differential equations: 

(6) = /(*)« + f{t)-u x + g{t)u t + g{t)u + xt 

x 

Taking the /^-transform of PDEs of the firm shown in ©, we get: 

- u(l + 2s s )f{t) + u t (l + 2s 2 )g(t) ^ = J(t) 

u 

Where J(t) = — ypy- This implies that a solution for u is as follows: 

u{ Sl t)=eti J{w)dw -L{s) 

Where L(s) is an arbitrary function of s. Representing this expression as a series, 
and letting L(s) = ^2, we get: 

^ n\ ■ 2y 2 

n=0 y 

And after applying the C 2 inverse to our sum, we get the following solution for 
u(x, t): 

u{x, t) = x° ■ eti J ( w ) dw ) 

5. The C 2 Convolution 

Definition 2. A binary operation fk) called the convolution of two functions f 
and g is defined as follows: 

(/*<?)(*) = [ xf(Vt 2 -x 2 )g(x)dx 
Jo 

It can be shown that this operation is associative and commutative. Most im- 
portantly for our purposes, the following are true: 

C 2 {f*g;s} = £ 2 {f)-£ 2 (g) 

and: 

(7) C 2 {h * h * ... * fn, s} = C 2 (fi) ■ C 2 (f 2 ) ■ ... • C 2 (f n ) 
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From which it follows that: 

(8) /THA ' h •••••/„} = fl *h *•••*/« 

6. An Application to a Partial Differential Equation that is of 
Exponential Squared Order, and not of Exponential Order 

Consider the following partial differential equation, with the following condition: 

(9) = g{t)u - f{t)u t + -f{x)u xt u{0+,t) = 

x 

Writing this equation in terms of the differential operator, 

= g(t)u - f{t)u t + 6 x f(t)u t 
and applying the /^-transform, we arrive at the following expression: 

= g{t)u t - f(t)u t + 2s 2 f(t)u t 

Rearranging, 

= g(t)u + f(t)f t (-l + 2s 2 ) 

Which leads to: 

Where H(t) — -|^y. This implies that a solution for u is as follows: 

u(s,t) = e-!+ 2 ^ Jo ^ ; 
Writing this solution as a series, we get: 

(j*H(w)dwr 



n 



n=0 

from ([7]), © and the following identity: 

1 



! \-l + 2s 



C 2 {e ax ;s} = 



-2a + 2s 2 

We can arrive at the following solution: 



n=l \ 

Where (e^ 1 / 2 ^ )*« represents the following: 



e (l/2)e^ e (l/2)e 2 _^ e (l/2)e 



Which leads to the following proposition, which can be shown through induction: 

„ ,„s 2 1 x An-i) . M/2)x 2 
( e (l/2)e 2 )*n = 1 . f e 



2"- 1 (n-iy. 

For n > 1. Applying this proposition to the cases in which n > 1, we have the 
following solution for u(a:,i): 

' (/„' H{w)dw) n 1 x 2 ^ 1 ) • e^ 1 / 2 )" 2 \ 



in — 1 



(n-1)! 
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We now make some remarks regarding the order of our solution. For the sake of 
simplicity, we assume that our function J Q H(w)dw — 1, giving us the following: 



(2n-2) . e - 



n=l 

Notice that: 



(x, i) = 1 + V — ^— 



« / x (2n-2) . g^ _ 

t(x, i) = 1 + > -— — = — t — = 1 + e~ ■ > 



, n! ■ 2™- 1 • (n - 1)1 / ' ^ (n + l)-(n!) 2 -2" : 

n=l \ X ' / n=0 y ' y ' 



X 2n - 



and that 

oo 

^Tn + V, ■ \n'r -2" 

n=0 V 1 / V / 

is a sum S of elements in R + , where one can take the sum of first two terms where 
n = 0, 1 as evidence that S > 1. From this it follows clearly that 

e~ ■ S > e~, 

and since 

lim e ^(1 + e 2 ) = oo, 
our solution is not of exponential order. Alternatively, notice that 
(n + 1) • (nl f ■ T > (n!) • 2" ■ 2", 

because ultimately, 

n\ > a n 

for some fixed a € R. From this inequality we conclude that: 



e 2 • > — — — — — — < e 2 • > — — — = e 2 • > - — '- — = e 2 • e 4 = e 4 

^ (n + 1) • (n!) 2 • 2" ~ ^ n\ ■ 4" ^ n 

n=0 v ; v ' n=0 n=0 



Noticing that 



lim e ^"(1 + e 4 ) = lim 



= lim e x + lim e 4 = 0, 



We see that our solution is in fact, of exponential squared order. Since our solution 
was not of exponential order, one would not have been able to calculate the solution 
to our PDE in equation ^ with a Fourier integral transform, which can only handle 
equations up to exponential order. 
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